Global persistence exponent of the two-dimensional Blume-Capel model 
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The global persistence exponent 6g is calculated for the two-dimensional Blume-Capel model 
following a quench to the critical point from both disordered states and such with small initial 
magnetizations. Estimates are obtained for the nonequilibrium critical dynamics on the critical line 
and at the tricritical point. Ising-like universality is observed along the critical line and a different 
value 6g — 1.080(4) is found at the tricritical point. 
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Recent advances in nonequilibrium physics have enlarged our understanding of critical phenomena. As shown by 
Janssen et al. [1] and Huse [2], universality and scaling behavior are already present in systems since the early stages of 
their relaxation processes after quenching from high temperatures to the critical one. This evolution is governed by a 
new independent exponent 9. Systems characterized by a nonconserved order parameter, model A in the terminology 
of Halpcrin et al. [3], are described by a scaling function soon after a microscopic time scale tmic- Its general form for 
the fcth moment of the magnetization (e.g., in a ferromagnet) reads as 



Here b is an arbitrary spatial scaling factor, t is the time evolution and r is the reduced temperature, t = [T — T^jT^. 
As usual, the exponents /3 and v are the well-known static exponents, whereas z is the dynamic one. Equation (1) 
depends on the initial magnetization mo and gives origin to the new exponent xq^ scaling dimension of the initial 
magnetization, related to by = 0z + iijv. 

The nonequilibrium short-time exponent can be obtained at the critical temperature from the scaling form of the 
first moment [1,4,5], M(t,mo) ~ mot^ . Usually, this exponent assumes a positive value corresponding to a critical 
initial slip, which is related to the anomalous behavior of the magnetization when the system is quenched to Tc. 
Numerical works indicate negative values for 9 at the tricritical point [6] of the S = 1 Blume-Capel model [7-9] 
in two dimensions {2d) as well as for the Baxter- Wu [10], Ising model with three-spin interactions [11] and for the 
four-state Potts model [12]. The indication of a negative value for this exponent was theoretically deduced by Janssen 
and Oerding [13] from a study of non-equilibrium relaxation at a tricritical point. Numerical simulations of the 
S = 1 Blume-Capel model [6] present for the dynamical exponents the values 6 = —0.53(2) and z = 2.21(2) at the 
tricritical point and values compatible with heat-bath dynamics for the two-dimensional Ising model {6 — 0.191(1) 
[4,14], z = 2.156(2) [15]) on the critical line. 

Under the same nonequilibrium conditions, a second critical exponent has been presented in the literature [16]: 
the global persistence exponent 9g. It is related to the probability P{t) that the global order parameter has not 
changed sign up to time t after a quench to Tc, P{t) ~ i^^f. This exponent has emerged from the concept of the local 
persistence phenomena in coarsening dynamics at T = [17-22]. As argued by Majumdar et al. [16], if the dynamics 
of the global order parameter is described by a Markovian process, 9g is not a new independent exponent and it can 
be related to other critical exponents, 



(1) 



0gZ = \-d+l-ri/2. 
Here A is the nonequilibrium exponent of the auto-correlation function [1,2], 
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which is related to the short-time exponent, X = d — 9z. Therefore, Eq. (2) can be rewritten as 



2 V 



(4) 



However, the time cvohition of the order parameter is, in general, a non-Markovian process and 6g turns out to be a 
new independent critical exponent describing the stochastic dynamic process toward the equilibrium. 

Contrary to the local persistence, the global persistence has been less studied. Results have only been reported for 
the n — !■ cx) limit of the 0{n) model {9g = (d — 2)/4 for 2 < d < 4 or 0g = 1/2 for d> A) and to order e = 4 — d near 
d = 4.{eg = l/2-e(n + 2)/(4n + 32)+0(e2)), for d = 1 Ising model {Og = 1/4) [16], and d = 2 Ising model [16,23,24], 
as shown in Table I. As remarked in [16], relation (2) is satisfied for n = oo limit of the 0{n) model, to first order 
in e = 4 — d, and also for d = 1 Ising model. However, for the 0{n) model, it has been shown [25] that the scaling 
relation (2) for Markovian process is violated at order e^. 

In this paper we investigate the universality aspects of the global persistence exponent for the Blume-Capel model. 
This exponent is obtained by the straight application of the power law behavior P{t) ^ t~^^ and by means of time 
series data collapse. This study is also performed under different initial conditions: random choices of all spins and 
sharp preparations of samples with defined and nonzero magnetizations mo [24] . 

The Blumci-Capel [7] (BC) model is a spin-1 model which has been used to describe the behavior of ''He — ^He 
mixtures along the A line and near the critical mixing point. Apart from its practical interest, the BC model has 
intrinsic interest since it is the simplest generalization of the Ising model (s = 1/2) exhibiting a rich phase diagram 
with first and second-order transition lines and a tricritical point. The Hamiltonian of the two-dimensional model is 



where < i,j > indicates nearest neighbors on lattices and Si = {—1,0,1}. The parameter J is the exchange 
coupling constant and D is the crystal field. We remark that along the critical line, this model presents a critical 
behavior similar to the Ising model. However, exactly at the tricritical point the exponents change abruptly. In [26] 
finite-size scaling combined with conformal invariance permitted to observe a smooth change between Ising-like and 
tricritical behavior. Ising-like behavior is reached only when L ^ oo, leading to the exact values of the Ising model 
critical exponents. 

In order to evaluate the persistence probability P{t), we define p{t) as the fraction of samples which change 
their signals for the first time at the instant t. Our probability P{t) is numerically calculated from the cumulative 
distribution function such that the total magnetization does not cross the origin up to time t, 



The spins {Si\ are updated by the heat-bath algorithm and our statistics rely on Nb = 5 independent bins with 
Ns = 40000 samples for t up to 1000 MC sweeps and lattice size L — 80. We quote estimates for the time intervals 
[ti,t/] with the highest values of goodness-of-fit Q [27] for the linear regression obtained at every 10th measurement. 

Table II lists the points on the second order critical line and the tricritical point where we have performed our 
simulations. This table presents estimates for 9g in function of different magnetizations mo to explore the effect of 
the initial configurations on the behavior of P{t). Linear extrapolations for mo — * are presented in the last column. 
Here we follow the sharp preparation technique to set the value mo. Our typical time intervals for mo = 0.0050 
correspond to [100,500] (critical points) and [40,400] (tricritical point). Different time intervals with accepted values 
for Q, present compatible results within our error estimates. 

Figures 1 and 2 illustrate, respectively, the decay of persistence probability for the specific critical point D/J = 0, 
ksT/J — 1.6950 and for the tricritical point. 

Our simulations of the BC model on the critical line reproduce (see Table II) the estimates obtained by Schulke 
and Zheng [24] for the 2d Ising model with mo = 0.0005: 9g = 0.238(3) (with the HB dynamics) and 9g = 0.236(3) 
(Metropolis algorithm). On the other hand, our simulations present the largest deviation (compared with the Ising 
model exponent) for the largest initial magnetization mo = 0.0050, but the expected universality is recovered as 
mo — > 0. It becomes clear the importance of the initial configurations in measuring the persistence exponent. 

At the tricritical point we observe a faster decay of the persistence probability {9g = 1.080(4)), characterizing a 
different persistence behavior. We repeat our simulations for two new critical points [28] D/J = 1.87, ksT/ J = 0.800 
and D/J = 1.95, ksT/J = 0.650, closer to the tricritical one to follow this cross over effect. In this case, the 
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persistence probability can be fitted by a power law function only for short time intervals. Hence, we have to restrict 
our analysis to shorter time intervals in order to obtain acceptable values for the goodness-of-fit. We obtain for the 
first point 9g = 0.172(3) in time interval [80,400], while for the second point, closer to the tricritical one, we have to 
restrict further the interval to [100,300], leading to the estimate 6g = 0.345(5). Hence, this very different behavior 
characterized by those new values of d indicates a 'transition' from the Ising-like 6g ~ 0.23 to the tricritical 6g ~ 1.08 
for finite lattices. 

The initial magnetization dependence of P{t) can be cast in the following finite-size scaling relation [16], 



P{t) = t-'^ f(t/L^) 



(7) 



which renders a different method to obtain the exponent 9g from lattice sizes Li and L2 [16]. For this end we define 
W{t,L) = L^!i^-P{t), which turns out to be a function of t/L^. Therefore, if we fix the dynamic exponent z, the 
exponent 9g can be obtained by collapsing the time series W{t2, L2) = f{t2/L2) onto W{ti, Li) = f{ti/Ll) as follows. 
Under re-scaling, with b = L2/L1, {L2 > Li), we obtain 



W{t2,L2) = W{bHubLi), 
and the best estimate for 6g corresponds to the minimization of 



(8) 



W{t,L) - W{bH,bL) 



\W{t,L)\+ W{bH,bL) 



(9) 



by interpolating W to the time values b^t. 

This method is applied to the critical point D/ J = 0, ksT/J = 1.6950 and to the tricritical one. We performed 
simulations with lattice sizes L = 10, 20, 40 and 80 and initial disordered states to study the finite size dependence 
of 9g. Our simulations also rely on Ng = 40000 samples and Nb = 5 bins. The collapse obtained from our largest 
pairs of lattices {Li,L2) = (40,80) is displayed in Fig. 3 for the tricritical point. A similar figure (not shown) is also 
obtained for the critical one. Results for the persistence exponent are presented in Table III with the input values 
z = 2.106 and z = 2.215, respectively for the critical and tricritical point [6]. Here, we can observe a good agreement 
between both methods. However, we remark that the collapse method gives less precise results. 

Now, we return to Eq. (4) to analyze the markovian aspects of the time evolution of the magnetization. We verify 
that our numerical precision for both the exponents [6] 9 and 9g, evaluated on the critical line and at the tricritical 
point, allows us to detect the non-markovian behavior. Thus, these results substantiate the independence of global 
persistence and short-time exponents, characteristic of a non-markovian dynamic evolution of the magnetization. 

In summary, we have studied the effects of the initial magnetization on P{t) and evaluated the global persistence 
exponent of the Blume-Capel model from different methods (linear extrapolation toq ~^ and collapse). The uni- 
versality and independence of this dynamic exponent is explicitly shown for the BC model along the critical line. 
In addition, the power law time dependence of the persistence probability is also exhibitted at the tricritical point, 
presenting a faster decay when compared with the Ising model exponent. 
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TABLE 1. The global persistence exponent 6g for the 2d Ising model with random and with pre-fixed small initial magneti- 
zation mo. 



Reference 


random 


mo = 0.0010 


mo = 0.0005 


[24] 


0.233(5) 


0.237(5) 


0.238(3) 


[23] 


0.225(10) 






[16] 


0.233(9) 







TABLE 2. The global persistence exponent dg from the power law behavior for the 2d Blume-Capel model for different initial 
magnetizations mo. Last column contains our linear extrapolation to mo = 0. 



D/J 


ksT/J 


mo = 0.0050 mo = 0.0025 


mo = 0.0005 extrapolated value 





1.6950 


0.219(3) 0.233(2) 


0.237(3) 0.241(4) 


-3 


2.0855 


0.223(2) 0.233(1) 


0.234(2) 0.238(4) 


-5 


2.1855 


0.223(2) 0.232(1) 


0.235(1) 0.237(2) 


tricritical point 








1.9655 


0.610 


1.054(3) 1.072(4) 


1.076(3) 1.080(4) 


TABLE 3. The ; 


global persistence exponent 6g for the best data collapse for the 2d Blume-Capel model with initial magne- 


tization mo = 0. 












critical point 


tricritical point 






{D/J = 0, fesT/J = 1.6950) 


{D/J = 1.9655, ksT/J = 0.610) 






0, 


O9 


10 


20 


0.254(7) 


0.88(2) 


20 


40 


0.26(2) 


1.08(2) 


40 


80 


0.24(2) 


1.06(2) 
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FIGURE CAPTIONS 



Figure 1: Persistence probability P{t) for the 2d Blume-Capel model for L = 80 at a critical point with the sharp 
preparation of the initial magnetization mo- 

Figure 2: Persistence probability P{t) for the 2d Blume-Capel model for L = 80 at the tricritical point with the 
sharp preparation of the initial magnetization mo- 

Figure 3: Collapse of persistence probability at the tricritical point D/J = 1.9655, ksT/J = 0.610. 
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FIG. 1. Persistence probability P{t) for the 2d Blume-Capel model for 1/ = 80 at a critical point with the sharp 
preparation of the initial magnetization mo- 
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FIG. 2. Persistence probability P{t) for the 2d Blume-Capel model for L = 80 at the tricritical point with the sharp 
preparation of the initial magnetization mo- 
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